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Abstract -In this paper we study the scaling behavior of the fluctuations in the steady state Ws 
with the system size N for a surface growth process given by the competition between the surface 
(y^ _ relaxation (SRM) and the Ballistic Deposition (BD) models on degree uncorrelated Scale Free 

^ ■ networks (SF), characterized by a degree distribution P(k) ~ fc _A , where k is the degree of a node. 

' It is known that the fluctuations of the SRM model above the critical dimension (d c = 2) scales 

logarithmically with N on euclidean lattices. However, Pastore y Piontti et. al. [A. L. Pastore y 
Piontti et. al, Phys. Rev. E 76, 046117 (2007)] found that the fluctuations of the SRM model 
C*j | in SF networks scale logarithmically with N for A < 3 and as a constant for A > 3. In this letter 

O f we found that for a pure ballistic deposition model on SF networks Ws scales as a power law with 

an exponent that depends on A. On the other hand when both processes are in competition, we 
find that there is a continuous crossover between a SRM behavior and a power law behavior due 
to the BD model that depends on the occurrence probability of each process and the system size. 
Interestingly, we find that a relaxation process contaminated by any small contribution of ballistic 
deposition will behave, for increasing system sizes, as a pure ballistic one. Our findings could be 
■^j- ■ relevant when surface relaxation mechanisms are used to synchronize processes that evolve on top 

t-H ' of complex networks. 

(N- 

■ In the last decade the study of complex networks received much attention because many 

k^j | processes take place on top of these kinds of structures. Here, we consider networks with 

scale free (SF) topologies characterized by a degree distribution Plk) ~ k~ x , where k is the 
degree or number of connections that a node can have with k rnax > k > k m i n , where k max 
and k m in are the maximal and minimal degree respectively, and A measures the heterogeneity 
of the distribution [TJ[2] . Historically, the research on complex networks was mainly focused 
on how the topology affects processes that evolve on top of them, such as epidemic spreading 
[3J, traffic flow IH[S], cascading failures [Hill] and synchronization jBrHU]. Recently, it was 
shown that the scaling behavior of processes, such as synchronization and jamming [PHlT^] . 
are fully governed by the topology of the underlying network. Synchronization processes are 
very important in many real situations such as supply-chain networks based on electronic 
transactions [T3] , neuronal networks [Tl] and financial transaction between traders [T^ ■ One 
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of the most successful attempts to model synchronization phenomenon is to map them into 
a non equilibrium surface growth |101ll6j where the relevant magnitude that represents the 
departure from synchronization is the dispersion or fluctuations of some scalar field h over the 



nodes of a network given by W — < l/N J2i=i(^i ~ C 1 ) ) 2 f > where hi represents the scalar 



on node i, (h) is the mean value, TV is the system size and {.} denotes an average over network 
configurations. Synchronization problems deal with the optimization of the fluctuations 
in the steady state, being the system optimally synchronized when those fluctuations are 
minimized. On the other hand a system is called scalable if their fluctuations do not depend, 
or depend weakly, on the system size. 

Many dynamic synchronization processes can be described by simple growth models. One 
of the most successful models of surface growth with relaxation is the Family model |17j 
(SRM) that represents very well film growth either by vapor or chemical deposition. Other 
very important model is the Ballistic Deposition (BD) [TB] that represents sedimentation 
processes such as low thin-film grow at low temperatures [T5]. However, more complex 
processes cannot be described by a single growth model and are better represented as the 
competition between two or more simple mechanisms. The competition between models on 
cuclidcan lattices [2TTH25] was studied by few researchers, in spite of the fact that it is more 
realistic in describing grow processes on real materials. Pellegrini et. al. |20j studied the 
competition between the SRM model with probability p and the BD model with probability 
1 — p on cuclidcan lattices and found a crossover between a logarithmic and a power law 
behavior of the fluctuations in the steady state Ws with TV that depends on p for dimensions 
d > 2. Note that d — 2 is the critical dimension d c of the SRM model above which the 
behavior of Ws takes the mean field value, with Ws ~ lnJV. Pastore y Piontti et. al. [53] 
studied the SRM model on degree uncorrelated SF networks and found that 



a very different scaling behavior than the one obtained on euclidean lattices above d Cl even- 
though complex networks can be thought as high dimensional systems with d 3> d c . These 
results were theoretically confirmed by La Rocca et. al. [25] who derived the evolution 
equation of this model on complex networks. However, up to our knowledge, neither the 
BD model nor the competition between different processes on complex networks was ever 
reported. Usually the SRM model is associated with diffusion between first neighbors of a 
network, such as in load balance problems in parallel processors, where the load are tasks 
to solve. If the load difference between processors is large, nodes with less load must wait 
for the more loaded nodes to finish their tasks. Thus it is logical to send all the new task 
to the less loaded node to prevent the desinchronization of the system. Then sending new 
tasks to the less loaded processors could be represented by the BD model. 

In this letter we study the scaling behavior of the fluctuations in the steady state Ws 
with the system size TV for the BD model, and the competition between the SRM and the 
BD models on degree uncorrelated SF networks. 

In our simulations, at each time step a node i is chosen with probability l/N. Then, the 
SRM rules are applied with probability p and the BD rules with probability 1 — p. Denoting 
by Vi the nearest neighbor nodes of i, the evolution rules for the SRM model are 





lnJV for A < 3 , 
const, for A > 3 ; 



(1) 



(1) if hi < hjij € Vi =4> hi = hi + 1, else ; 

(2) if hj < h n in ^ j £ Vi hj = hj + 1 . 



Meanwhile for the BD model the evolution rules are 



(1) if hi > hjij e Vi hi = hi + 1, else , 

(2) if h n > hjij 7^ n € Vi => hi = h n . 
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As initial condition we choose {hi} = for i = 1, N and the time step increases by 1/N . 
The SF networks with different values of A were constructed using the configurational model 
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Fig. 1: Fluctuations Ws as function of N in log-log scale for pure BD model (p = 0) in SF 
networks for A = 2.5 (O), 2.8 (□), 3.0 (o), 3.2 (A) and 3.5 (*). The dashed lines correspond 
to a power law fitting for large system sizes N with Ws ~ N a . 



In Fig. [TJwe show the fluctuations Ws as a function of N for the BD model (p = 0) in 
log-log scale for different values of A. We can see that for all the values of A there exists a 
power law behavior denoted by the dashed lines that represent the fitting of the data with 
Ws ~ N a for large system sizes, where the exponent a depends on A. Clearly, we can 
observe that this behavior is very different than the one found for the SRM model on SF 
networks (see Eq. (TTJ). The power law behavior of the fluctuations in the steady state of 
the BD model seems to be an inherent property of the model regardless of the topology of 
the complex network over which the process evolves, although the exponent of the power 
law depends on the heterogeneity of the network. We find that as the heterogeneity of the 
network decreases, above A = 2.5, finite size effects become more relevant modifying the 
power law behavior found for large N. Thus, we propose that 

Ws ~ N a (l + AN~ S ) , (2) 

where the term N~ 5 is the correction to finite size scaling with S > and A a constant. We 
can measure the finite size effects taking the derivative of Ws /N a with respect to N 

x(A0 = A(^V^- (6+1) - 

XK ' 8N\N a J 

In Fig. [5] a) we plot x(iV) as a function of N up to N = 25000, for two values of A > 2.5. 
The dashed lines represent the fitting with a power law from where we obtain the values 
of 5. In order to show that our proposed scaling (see Eq. ©) is correct in Fig. [3]b) we 
plot Ws/N a as function of N~ s . We can see that for large system sizes Ws/N a does not 
depend on N , but as N decreases we observe a linear behavior, which shows the agreement 
between the proposed corrections to scaling and the simulations. 

Next, we study the competition between the SRM and the BD models, i.e., the scaling 
behavior of Ws with N when we increase the parameter p. Although the conclusions reached 
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Fig. 2: a) Corrections to the scaling x(N) as function of N in log-log scale for A = 2.8 (O) 
and 3.5 (□). The dashed lines correspond to the power law fitting from we obtain 5. b) 
Linear-linear plot of Ws/N a as function of N~ 5 for the same values of A than in a). 



from the competition between the two models are qualitatively the same for all the values 
of A, we will consider only the case A = 2.5 because the effects found are more pronounced 
in this case and because for this value of A finite size effects are negligible as we mention 
above. In Fig. [3] we plot Ws as function of N in log-log scale for different values of p. Wc 
can see that for p = there is a pure power law with exponent 0.186 (see FigQ} and for 
p = 1 the scaling behavior of W s with N corresponds to a logarithm as expected for A < 3 
[24l[25]. For values < p < 1 we find that there is a crossover between both regimes. The 
corresponding curves become steeper as p increases leading to expect that they will cross at 
some point. However, as the system size is increased we observe that the curves change their 
slope, avoiding the crossing between them. This is can be expected since the fluctuations in 
the competition with the SRM model, which tend to smooth out the interface via diffusion, 
cannot overcome the fluctuations of the pure BD model. In fact, we expect that for large 
sizes (N — > oo), the behavior of all the curves with p < 1 will tend to follow a power law 
with the same exponent that for the case p = 0. However from our data we cannot show our 
expectation even for the system sizes as big as the ones that we have achieved (N ~ 10 6 ). 
Note that each point in Figs. Q] and [3] corresponds to 10 4 realizations of the networks for 
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fix N where each value of Ws was obtained by a linear fitting on the steady state. Thus 
to reach bigger system sizes is very time consuming and goes beyond our computational 
capabilities. For this reason from the results shown until here we can only conclude that 
in the competition between both models, even for p close to one, the BD dominates the 
behavior of the fluctuations. 




Fig. 3: Fluctuations Ws as function of N in log-log scale for SF networks with A = 2.5, 
kmin = 2, for p = (O), 0.6 (□), 0.8 (o), 0.9 (A), 0.95 ( V ), and 1 (★). 
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Fig. 4: Schematic of the formation of steps for a one-dimensional lattice where we use the 
threshold Ah = 2. The steps, represented by different colors and capital letters, are bounded 
by nodes whose height difference exceeds the value of 2. The numbering in the horizontal 
axis represents the node number and the vertical axis is the height of each node. 



In the competition, the BD model is responsible of generating bigger height differences 
(jumps) between neighboring nodes than the SRM model, which tends to smooth out the 
interface by diffusion (see Fig. [3]). As a consequence, the big jumps that produce the 
BD model increase the fluctuations of the system. As N increase, the number of jumps also 
increases and as a consequence the BD model dominates the behavior of the system, even for 
p close to one. In order to verify our assumption we introduce the definition of "step" . Steps 
are landing structures with edges generated by the BD process and flat regions smoothed 
by the SRM process. Since Ws represents the dispersion of the heights around their mean 
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value, we can assume that most of the nodes of the network have a height /jk (/i)± mWs, 
where m is a positive integer number that determines the degree of confidence for a Gaussian 
distribution, i.e., m = 1 correspond to 66% of confidence and m = 3 to 99%. We verified that 
the distribution of heights is a Gaussian distribution for any value of p. This implies that in 
average the absolute value of the difference of heights between any two pair of neighboring 
nodes is at most 

Ah « V2mW s ■ 

Then, if we denote the fluctuations for the pure SRM model (p = 1) as W§ RM , we define 
a step as an almost flat connected region of the interface where each pair of neighboring 
nodes has Ah < \/2mW§ RM . Thus, computing how many steps each model generates and 
how this magnitude changes with p, we can measure the influence of the BD process [37j . 
In Fig. @]we show for Ah > 2 a schematic plot for the particular case of a one-dimensional 
lattice where each step is marked with a different color and labeled by a letter. For example, 
the step C consist of nodes IV and V that have Ah < 2 between them, and is limited by 
nodes III and VI because their height difference is greater than 2. Note that a step can 
be formed by a single node, as in the cases A and D. Thus a step is a soft region of the 
interface bounded by two borders represented by nodes with an abrupt jump in their height 
differences (Ah > y/2mW§ RM ). 
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Fig. 5: Average number of steps n$ as a function of N in the surface for m — 3 in log-log 
scale for SF networks with A = 2.5 and p = (O), 0.6 (□), 0.9 (o) and 1 (A). The lines are 
used as guides for the eyes. 

In Fig. [5] we plot in log-log scale the average number of steps ns as a function of N 
for p = 0, 0.6, 0.9 and I, for m = 3 and taking W§ RM 1 (see Fig. [3]). We can see that 
for p = 1, ns is much smaller than for p = for all the values of N. For large N and 
< p < 1 the curves of ns try to reach the value of p = 0, which means that the BD is the 
dominating mechanism (see Fig. [3]). Thus even for values of p very close to 1, corresponding 
to pure SRM, as N increases more steps appear and the BD model ends dominating the 
behavior of the fluctuations in the steady state. This explain the crossover behavior of Ws 
with N for all the values of p < 1. Thus a relaxation process driven by diffusion competing 
with a small contribution of the ballistic deposition will behave as a pure ballistic one for 
large N. Then, if we are interested in reducing the fluctuations in order to improve the 
synchronization of the system when these two processes are competing, we must increase 
the diffusion and reduce the system size. Our findings could have a great impact when 
surface relaxation mechanisms are used to synchronize system on complex networks because 
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a small contribution of the BD process reduces the synchronization and as a consequence 
the scalability of the system with the system size worsens [28] . 

In summary, we studied the scaling behavior of the fluctuations on the steady state for 
the ballistic deposition model in degree uncorrelated SF networks and found that Ws ~ 
TV" where the exponent a increases with the heterogeneity of the SF networks, i.e., as A 
decreases. Our results suggest that the power law behavior of Ws is an inherent property of 
this process found also in euclidean lattices. We also studied the scaling of Ws with N for the 
competition between the SRM and the BD models and found that there exists a crossover 
between different behaviors as we increase p. For large N, we found that for p < 1 the 
values of Ws approaches to the BD one, which means that the BD model always dominates 
the behavior of the fluctuations in the steady state. This was shown using the definition of 
step, where the average number of steps measure the influence of the BD process. We found 
that the crossover between the SRM and the BD corresponds to an abrupt increase in the 
average number of steps. 

In a forthcoming work we will study analytically the BD model and the competition 
between the BD and the SRM models deriving the evolution equations for the heights from 
where we will be able to derive the scaling behavior of Ws for all the values of p. 

* * * 
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